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Abstract 

We study the excitation energy spectrum in the 5 = 1/2 ferro- 
magnetic Ising spin chain with the easy axis z in a magnetic field 
h = {h x ,0,h z }. According to Wu and McCoy's scenario of weak 
confinement, the fermionic spinon excitations (kinks), being free at 
h g = in the ordered phase, are coupled into bosonic bound states at 
arbitrary small h z > 0. We calculate the energy spectrum of such ex- 
citations in the leading order in small h z , using different perturbative 
methods developed for the similar problem in the Ising field theory. 

1 Introduction 

In recent years much progress has been achieved in the understanding of 
various aspects of the scaling limit of the two-dimensional Ising model, which 
is known as the Ising Field Theory (IFT), for a review see [4]. It is described 
by the Euclidean action 

•Aift = -Acft + t J e{x)d 2 x + h J a{x)d 2 x. (1) 



1 



Here Aift corresponds to the conformal field theory in two dimensions, which 
is associated with the the Ising critical point. Fields e(x) (energy density) 
and cr(x) (spin density) are the only relevant operators in the theory, their 
scaling dimensions are X e — 1 and X a = 1/8. Parameters r and h are 
proportional to the deviation of the temperature T and magnetic field H in 
the lattice Ising model from their critical values: r ~ T c — T and h ~ H at 
T — > T c , H — > 0. In fact, only one dimensionless parameter 77 = r/|/i| 8//15 
determines the physics of IFT. 

IFT is integrable along the directions h = and r = 0. The former case 
h = 0, r 7^ corresponds to Onsager's solution [16]. Exact solution of IFT 
at r = 0, h 7^ was found by A.B. Zamolodchikov [23]. A thorough study 
of analytical properties of the free energy continued to complex values of the 
scaling parameter 77 has been done by Fonseca and A.B. Zamolodchikov [10] 
by means of a numerical technique, known as the truncated free-fermion 
space approach [22]. Their analysis clarifies the relationship between the 
edge Yang-Lee singularity and the spinodal point. 

For an analytical study of IFT for h and r close to the integrable di- 
rections, it is natural to exploit perturbation expansions. Form-factor per- 
turbation theory developed by Delfino, Mussardo and Simonetti [6] has been 
applied to calculate the variation of the particle mass spectrum and the decay 
widths of non-stable particle for small 77, i.e. near the line r = [5,6]. 

Perturbation theory for IFT in the region close to the h = axis turns out 
to be rather nontrivial in the low temperature phase T < T c , as it was first 
shown by McCoy and Wu [14]. At h — the particle sector of IFT contains 
one spinless fermion, which is interpreted in the ordered T < T c phase as a 
kink interpolating between two degenerate vacua. Application of magnetic 
field removes the degeneracy and induces a long-range interaction between 
kinks, which leads to their confinement into pairs. This means, that even 
at small field h an isolated kink gains an infinite energy, and the bounded 
kink-antikink pairs become the only single-particle excitations in the model. 
Their dispersion law E n (P) has the relativistic form: 

E n (P) = y/P* + M*, (2) 

being completely determined by the masses M n . This follows from the 
Lorentz invariance of IFT, which in turn is the result of the rotation in- 
variance of the two-dimensional Ising model in the scaling limit. 

At small h the mass spectrum M n of kink-antikink pairs becomes dense 
in the segment [2m, 00). Two asymptotic expansions describe M n at h — > 
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in different regions of this segment. Near the edge point 2m (i.e. for fixed n 
at h — > 0) one can use the low energy expansion in fractional powers of the 
magnetic field 

i§ = !+f>C" 3 , (3) 

k=2 

where m = 2nr is the mass of a free fermion (kink) at h — 0, a is the spon- 
taneous magnetization at h = 0, and ( = 2ha/m 2 . Several initial coefficients 
fik are known explicitly. The coefficient of the leading term has been found 
by McCoy and Wu [14] in the form /i 2 = z n , where — z n denotes the zeroes 
of the Airy function, Ai(— z n ) = 0. Correction terms in ([3]) up to k = 8 have 
been calculated by Fonseca and A.B. Zamolodchikov, see [8,10]. 

On the other hand, for n 3> 1 and h — > the semiclassical expansion in 
integer powers of h can be applied: 

M 2 

_iL = [l + a 2 C 2 + 0(C 3 )]coshtf n , (4) 

where [9] a 2 = 0.071010809 . . ., and the numbers denote the solutions of 
equation 

sinh2tf n - 2$ n = 2ttC (n - 1/4) - C 2 5iK) - 0(C 3 ), (5) 

with 
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(6) 



The leading terms of order ( in (jSJ) were determined by Rutkevich [18], 
and independently by Fonseca and A.B. Zamolodchikov [8]. The second-order 
term ([6]) in dSJ) has been found by Fonseca and A.B. Zamolodchikov [8]. 

IFT viewed as a particle theory gives us a nice model of quark confine- 
ment. It is worth to note, that there are a lot of similarities between con- 
finement in IFT and in 't Hooft's model for the two-dimensional multicolor 
QCD [11], see the discussion in [8]. On the other hand, IFT can be treated as 
the continuous model of the one-dimensional quantum ferro- (or antiferro-) 
magnet. The continuous approximation applies to one- dimensional ferro mag- 
nets only in the critical region near the quantum phase transition point [20]. 
Beyond the critical region the discreteness of the lattice become important. 
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The purpose of the present paper is to discuss, how the picture of con- 
finement outlined above is modified by the discrete lattice effects beyond 
the critical region. Preliminary results in this directions have been obtained 
in [19]. Our work is motivated by recent experimental observation of confine- 
ment of topological excitations in spin-1/2 chain antiferromagnet reported by 
Kenzelmann et al. [13]. Accordingly, we shall hold to the magnetic terminol- 
ogy: the terms "spinon", "domain wall" or "kink" will be used for topological 
excitations, and their bound-states will be called "bound-spinons" [13]. The 
problem will be considered within the simplest appropriate one-dimensional 
model of ferromagnetism: S — 1/2 Ising spin chain in a "skew" magnetic 
field [3,7], which has components both normal and parallel to the easy mag- 
netization axis z, h = {h x , 0, h z }. This model has an exact solution in the 
case of a purely transverse magnetic field h = {h x , 0, 0} and reduces to IFT 
in the critical region near the quantum phase transition point [20]. We focus 
our attention on the dispersion law of bound-spinons in the limit of small 
longitudinal magnetic field h z — * 0, and extend the results (j3J), (jH) to this 
model beyond the critical region^. It turns out, that the bound-spinon energy 
spectrum in the discrete spin chain exhibits a rather rich small-/i z asymp- 
totical behavior. We describe eight different regimes for the bound-spinon 
energy spectrum, which are realized in the noncritical Ising spin chain in the 
limit h z — > 0, instead of two regimes ([3]), (jlj) known in IFT. 

The rest of the paper is organized as follows. In Section [2] the spin- 
1/2 Ising chain model in a skew magnetic field is described. We study the 
bound-spinon spectrum in this model by means of two different procedures, 
which are shown to be effective in the IFT [8,18]. In Section [3] we apply 
a heuristic approach based on the Bohr-Sommerfeld quantization rule. Re- 
sults obtained in Section [3] are confirmed and extended in more systematic 
theory in subsequent sections. In Section H] we derive a singular integral 
equation, which determines the energy spectrum of bound-spinons in the 
two-fermion approximation, and generalizes to the discrete spin chain model 
the IFT Bethe-Salpeter equation [8,10]. Section [5] contains some technical 
results relating to the singular integral equation (T4"3"l) . In Section Owe derive 
exact solutions of two simplified versions of this equation, and give explicit 
expressions for the corresponding bound-spinon energy spectra. In Section 
[7] we return to the original integral equation fj4*3|) . and develop for it a weak 

1 Of course, the relativistic functional form @ does not hold for the bound-spinon 
spectrum in the discrete spin chain model apart from the critical region. 



4 



coupling perturbative expansion based on the results of Section The lowest 
order of this expansion leads to eight different small- h z asymptotics for the 
bound-spinon dispersion law, which can be used in different regions of the 
bound spinon quasimomentum and energy. Obtained results are discussed 
in Section [HI 



2 Ising model in skew magnetic field 

The quantum Ising spin- 1/2 chain model is defined by the Hamiltonian: 

U=- J>] <J z j+1 + h x <r* + h z af). (7) 

3 

Here a x ' z are the Pauli matrices, j enumerates the chain sites. The number of 
the chain sites is put to infinity in the thermodynamical limit. The applied 
"skew" magnetic field h has both components parallel and normal to the 
magnetic easy axis z. The Hamiltonian is normalized to the ferromagnetic 
nearest neighbor coupling constant. 




critical region 



para 



ferro 
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Figure 1: Phase diagram of the Ising spin chain at T = 0. 



The phase diagram of model (j7j) at zero temperature is shown in Fig. [TJ 
Model (j7j) is not integrable at generic nonzero h x , h z . However, in the line 
h z = it has exact solution first found by Pikin and Tsukernik [17]. At h x = 
±1 lie the quantum phase transition points, which separate ferromagnetic 
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and paramagnetic phases [20]. In the critical regions near these points the 
model is equivalent to IFT. Two ferromagnetic ground states | 0+) and | 0_) 
coexist in the ferromagnetic phase — 1 < h x < 1. These ground states have 
the same energies E + {h z = 0) = E_(h z = 0), but the opposite signs of the 
spontaneous magnetization (0± | ex? | 0±) = ± a, where a — (1 — h 2 x ) l l s . 

In the integrable line h z = the model Hamiltonian ([7]) can be reduced 
to free fermions by use of the Jordan- Wigner transformation [12, 17]: 

— oo(6)a'(6)a(6) + const, (8) 

where 6 is the quasi-momentum, fermionic operators a^(6), a(6) satisfy the 
canonical anticommutation relations 

{a(6) ,a(9')} = {a f (0) ,a\6')} = 0, 
{a\Q) ,a{6')} = 2*5(6-6'), 

and the free-fermion dispersion law is given by 



u,(0) 



h x ) 2 + 4h x sin 2 - 



1/2 



(9) 



Operators crj and crfcf+i are bilinear in fermionic fields. Operator cr| can 
be written as a normally ordered exponential of a bilinear form of fermionic 
operators [12]. It can be also completely characterized by the form-factors: 

(9 1 , . . . , 6 K \a z \6[ ...,6' N ) = (OKft), . . . , a(6 K ) a z a\6[), . . . , a+(^)|0> 3 

where (K + N) takes even values. The elementary form-factors are given by: 

1 - exp^e*! + 6 2 )\ ^00(61)00(62) 

exp[-«(6li + 6 2 )\ - 1 ^00(6^00(62) 

a-\6\a^\6') = F(6\6>) = i-_ + (12) 

S ' 01 ' 1 ' ^ l-exp[i(0-0')] y/u{0)u{0') 1 ^ 

All other form-factors can be obtained from the elementary ones by applica- 
tion of the Wick rule. For example, the form-factor with two fermions in the 
initial and final states can be written as: 

a 
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(13) 



3 Bound states of two domain walls 

In the ferromagnetic phase \h x \ < 1, the free fermions in (J7j) represent 
the domain walls, which separate regions with different orientations of the 
magnetization. A small longitudinal magnetic field h z > 0, which breaks 
the Z2-symmetry, evidently provides the long range attraction force x = 
2h z a between the two neighboring domain walls. Due to this confining 
force, all domain walls become coupled into pairs which we shall call bound- 
spinons, following [13]. If h z is small, the weak confinement regime is realized. 
The energy spectrum of bound-spinons in this regime can be understood 
to much extent in the following heuristic approach, first developed for the 



Consider two interacting fermions moving in a line as a classical system 
with the Hamiltonian 



Here the fermion coordinates x±, %2 are the real numbers. Variables 6\, 62 
are the canonical momenta corresponding to the coordinates X\, X2, and the 
fermion kinetic energy uj{9) is given by After the canonical transforma- 
tion 



IFT [8, 10, 14, 18] . 



H(9%, 8 2 , x 1} x 2 ) = u(0i) + uj{9 2 ) + x \%2 - %i 



(14) 



X 



Xi + x 2 



2 



X — X2 — 



e 



#1 + 2 , 



02-e 
2 



(15) 



the Hamiltonian (j!4p takes the form 



H(8,x;Q) = e{6]Q)+x\x 



(16) 



where 



e(6; 9) 



cu{e + Q/2) + u(9-Q/2). 



(17) 



The total energy-momentum conservation laws read as: 



e(0(t); 9) + x \x{t)\ = E = Const, 
9(t) = Const. 



(18) 
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Figure 2: Function e(9] 9) given by (ITT)) . ([9]) for = and = 0.5. For 
energy E classically allowed region is [— 9 a ,9 a ]. 



The canonical equations of motion are: 



0(t) =-x signet)]. (21) 



For a given value of the conserved total momentum 0, Hamiltonian ( fT6l) 
describes the relative motion of two fermions. This motion becomes especially 
simple in the 6*-space. 

Fig. [2]shoes the ^-dependence of the "kinetic energy" e(9] 0) for h x = 0.5 
and = 0. If 9 — 9 a and x = +0 at t — 0, the total energy takes the 
value E = e(8 a , 0). Due to (12T1) . the momentum 9 of the relative motion will 
linearly decreases in time 

9{t)=6 a - x t (22) 

until the moment t — t\ — 28 a /x, when 9(t) reaches the value — 9 a . The time 
variation of the space coordinate x(t) for < t < t\ can be read from (TT8T) : 

= E- A e(t)-M 

X 
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Then x(t) becomes negative, and in the interval ti < t < 2ti canonical 
coordinates vary as 



9(t) = -9« + (t-ti)x, 
E- mt );B] 

X 

Then the phase trajectories periodically repeat in time with the period 2t\. 
The semiclassical energy levels can be obtained from the Bohr-Sommerfeld 
quantization condition: 

j dx9 = 2tt(z/+ 1/2), (24) 

where integration is taken along the closed periodic phase path in the (x, 8)- 
plane over one cycle of motion, < t < 2t\. Since the two interacting 
particles are fermions, only odd values are allowed for the integer v due to 
the Pauli principle: 

v = 2n-l, n= 1,2,... . (25) 

The right-hand side of (1241) can be transformed as 

dx6 = -2 [ a d99^P- = - j <W0:r(t)= (20) 
2 f\ e9 ^0;e) 




x J-». ae 

Here in the second equality we have taken into account the linear dependence 
([22]) between 6 and t for < t < ii, and in the third equality we have 
used equation of motion ([20]) . Combining (T2~4"|) . (1251) . and (1261) we find the 
semiclassical energy spectrum E n (Q) of bound states in the system of two 
fermions (fT4|) : 

2E n (Q)8 a - d0e(0,e) = 27rx(n-l/4), n=l,2,..., (27) 

where 6 a G [0, 7r] is the solution of equation 

e(6 a] Q) = E n (Q). (28) 
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In the above treatment we implied that function e(9; 0) monotonically 
increases with 9 in the interval < 9 < ir, providing that equation (1281) has a 
single solution there for E £ [e(0; 0), e(7r; ©)]. This is true for small enough 
values of the total momentum 0: 

< < m , m = 2 arccos /i-E. 

However, for > m the second derivative d 2 e(9; Q)/d9 2 at 9 = becomes 
negative, see Fig. [31 Though for high enough energies e(0, 0) < E < e(7r, 0) 
we can still use formula (1271) . for smaller energies E < e(0, 0), equation 
([15]) has two solutions 9 a , 9b in the interval < 9 < n. In the latter case 
the classical phase trajectory 9{t) moves back and forth in one of the the 
classically allowed region, say in the interval [9b, 9 a }. Equations (|2"2"j) . (|2"31) 
remain still valid, with t\ = (9 a — 9b) /x, however. Semiclassical energy levels 
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Figure 3: Function e(9; 0) determined by ( 1T71) (EI) for = 0.957T and /i^ = 
0.5. For energy E classically allowed regions are [— 9 a , —9b] and [9b, 9 a }. 

are determined by condition (|2^1) . where the integration path is located now 
in the interval [9b, 9 a }. Since this interval is not invariant under reflection 
9 — > —9, the Pauli principle does not impose any restrictions on the integers 
v = 0, 1, 2, ... in ([23). Therefore, instead of f[2T|) . we get for m < < tt 
and £ < e(O,0): 

E n (Q)(9 a -9 b )- d0e(0,e)=7rx(n-l/2), n = l,2,..., (29) 
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where n = v + 1, and 9 a , 9b E [0, 7r] are the solutions of equation 

6(00,6,9) = £7„(e), (30) 
shown in Fig. [3l 

Semiclassical formulas ff271) . (T291 can be used, if n is large, n ^> 1, and 
the energy E n (Q) does not approach one of the critical values E C (Q) of the 
function e(0,9): E e (Q) = e(0 c ,6), de(8,G)/d6\ e=Bc = 0. On the other 
hand, the low-energy part of the spectrum for small h z — > and < < O m 
can be determined by use of the "nonrelativistic approximation" developed 
for the IFT [8,10,14]. One should expand the "kinetic energy" in 9 near the 
origin e(9, O) ~ £o(9) + ^2(©)0 2 /2 reducing the problem to the Schrodinger 
equation with the Hamiltonian £o(9) — ^82(0)8^ + x\ x \- Only odd wave 
eigenfunctions <j>(x) should be taken into account due to the Pauli principle, 
and the energy levels take the form 

E n (Q) ^e (e) + x 2/3 h(O)/2] 1/3 ^, (31) 

where 

e (9) = e(0, O), e 2 (6) = d 2 e(9, Q)/d9% =0 , (32) 
and —z n , n = 1, 2, . . . are the zeros of the Airy function Ai(— z n ) = 0. 



4 Two-fermion approximation in the Ising chain 
model 

The bound-spinon energy spectrum can be formally defined as the solution 
of the eigenvalue problem 

(H - £ vac ) I $„(©)> = E n (Q) I $„(©)), (33) 

where 7i is the Ising chain Hamiltonian ([7]), E vac is the ground state en- 
ergy, and |$ n (0)) is the bound-spinon eigenvector with quasimomentum O: 
Ti|$ n (0)) = exp(iO) I $„(©)), with T x being the lattice site translation oper- 
ator. In the small-/i 2 limit this problem can be studied within the two-fermion 
approximation, which turns out to be effective in IFT [8,10,18]. In this ap- 
proach one considers instead of the the exact eigenvalue problem (|3~T!) its 
projection to the two-fermion subspace JF 2 of the Fock space 

V 2 UV 2 I $ n (9)) = £ n (9)P 2 I $ n (9)), (34) 
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where P 2 is the orthogonal projector onto JF 2 . So, the bound-spinon state 
is approximated by a two-fermion state, neglecting the four-fermion, six 
fermion, and higher multi-fermion contributions. In the momentum rep- 
resentation flMl) reads as 



[w(0O + u>(9 x ) - E n (Q)]U0u 82) = y jf ^^M8'i, d' 2 ) ■ (35) 

Wi + ®2 + - 0) (0 2) B l \a z \9[9 2 ), 

k 

where 



h,9 2 ) = (9 2 ,9 l \V 2 \<f> n (e)), 

0n(02,0l) = -0n(01,0 2 ), (36) 

0! + 2 = 6 + 27rfc, A; = 0, ±1, ±2 

Taking into account (1131) and (I36p , we can replace the two-fermion form-factor 
in the integrand, as 

(0 2 ,0i|^|0i0 2 )^4a^(0 2 ,0 1 |0i,0 2 ), 
^(0 2 ,0 1 |0i,0 2 ) = i[F(0 2 ,0 1 |)F(|0' 1 ,0 2 ) + 2F(0 1 |0' 1 )F(0 2 |0 2 )], 

where the elementary form-factors F(. . .) are given by fflUj) . Since the bound- 
spinon quasimomentum B is fixed, one can rewrite ( l35i) in the variables 

e = e 2 -e 1 and 0' = 2 - e[ ■. 

—G e (0, 000^0(0'), (37) 
where e(0; 9) is given by (|T7|) . x = 2/i z cr, and 



G (0,0') = S( | + ^| -0 



§ - f, | + f, ) (38) 



The kernel has the second order pole singularity at = 0': 

GM * = ~ {e-p w - y)/a| -eM-«e - + G ^ W ^ (39) 
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where G 



c> eg \ 



, $') is regular at real 9 and 9'. 
The function (p n ,e(9) in equation (1371) is odd and 27r-periodic in 9, and 
the integral in the right-hand side is understood in the sense of the Cauchy 
principal value. Equation (13"T|) gives the discrete-lattice version of the IFT 
Bethe-Salpeter equation in a generic momentum frame, see (3.11) in [8]. It 
is instructive to rewrite (l37j) in the coordinate representation: 



(40) 



K(j - f) n , e (j") - E n {Q) 4>n,eU) 



-X \j\ 0n,e(j) + X $^W(j, j') (f> n ,e{j'), 
f 

where K(j) and <p n ,&{j) are the Fourier coefficients of the 27r-periodic func- 
tions e{6; 9) and (p n ,e{9), respectively, = 0, ±1, . . .. 

The first and the second terms in the right-hand side of fj4*0|) come, respec- 
tively, from the first and the second terms in fl39|) . The first term describes 
the long-range attraction between fermions (domain walls), while the second, 
being responsible for the short-range interaction, vanishes exponentially in 
\j\ and on the correlation length scale. 

In this paper we shall focus our attention on the perturbative solution of 
the singular integral equation (1371) in the leading order in the small parameter 
X- As we know from the IFT, the short-range interaction does not effect the 
fermion bound-state energy in the leading order in h z , contributing only to 
the higher order corrections [8,18]. Accordingly, in what follows we shall drop 
the second term in the right-hand side of fl39|) . which does not contribute to 
the bound-spinon energy E n (Q) in the leading order in x- Then, equation 
f l37|) takes the form in variables z = exp(i9), z' = exp («#'): 



[e(z;Q)-E{e)]<K*) = -X'-t 



Si 



m (z — z) 



(41) 




hr+- ZV 



ZV 



1/2 



1 Z V 

h x + - 

ll x V Z 



1/2- 



(42) 



and v = exp(z0/2). We use notation S r for the circle of the radius r centered 
in the origin and passed in the counter-clockwise direction, and skip indices 
at E n (Q) and <j) n> e{z). 
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5 Singular integral equation in the unit circle 

In the present Section we describe some general properties of equation (I4ip . 
which will be used in Section [6] and UJ 

Consider a more general singular integral equation^] 



Functions e(z) and (f>(z) are supposed to be analytical in a narrow ring T, 
T = {z £ V \1 — 5 < \z\ < 1 + 5}, with small positive 5. The symmetry 
properties e(l/z) = e(z) and <j)(l/z) = —(f>(z) are also required. 
Let us define two functions g+(z) and g~(z): 



where g+{z) is defined at \z\ < 1, and g~{z) is defined in the region \z\ > 1. 
The evident properties of these function are: 

1. g+(z) and g~(z) are analytical at \z\ < 1 and at \z\ > 1, respectively. 

2. g+(z) and g~(z) can be continued to the unit circle Si, where they are 
continuous together with their derivatives. 



3- g+O-fz) = g-(z). 

4. s + (0)=s_(oo) = 0. 

5- g+(z) - g-(z) = (j){z) for \z\ = 1. 

6. d z9+ (z) + d,g.(z) = (vrz)- 1 ^ dz' <f>(z') (z' - z)~ 2 for |^| = 1. 
Combining properties [6] with fj43l) . one finds for \z\ — 1 

{[e(z) - A] +fizd s }g + {z) = U{z) = {[e{z) - A] - »zd z }g-(z). (45) 
The function U (z) defined this way is invariant with respect to the inversion 
U(l/z) = U(z), (46) 




(43) 




(44) 



2 On the general theory of singular integral equations see [15]. 
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and can be analytically continued from the circle Si into its neighborhood V. 
U(z) can be expressed directly in terms of function (f)(z): 

Note, that the integrand in (1471) is regular at z' — > 2. 

Let us prove this useful relation, following ideas of Vekua [21]. We start 
from the Cauchy formula 

(«) 

J ar 2m t- z 

where |z| = 1, and <9r = S14.5 — Si-a, ^±5 denote circles of radius 1 ± 5 
passed counter clockwise. Substitution of ( j4"5l) into ( )48j) yields 

crto = I — W)-x\9-{t)-M-{t) _ m 

dt [e{t)-X]g + (t)+fj,tg' + (t) 
Sl _ s 2ni t-z 

Several terms vanish in (j4"9l after integration, since g~(t) is analytical at 
\t\ > 1, and g+(t) is analytical at \t\ < 1. After substitution of (jH| in fj4"9]) . 
the result reads as 

7 ar 2m t- z J Sl 2iri z' -t 

Change of the order of integration leads to (J47J), since 

r dt e{t) e{z) - e(z') 

J av 2ni (t-z)(z'-t) ~ z'-z 

for z, z' G 5*1. 

It is helpful to inverse relation fj4"5l) and to express g+(z) in terms of U(z): 

g + (z) = f — t U(z') exp { - [F(z') - F(z))\, (50) 
J z t LZ [A 4 J 

where 

r dt 

iT{z) = J j[e(t)-X]. (51) 
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The choice of the integration path and its initial point Zq in (1501) depends on 
e(z), as will be discussed later. The general relations described here will be 
used in subsequent sections in analysis of equation (J4"3"j) for particular choices 
of function e(z). In Section [6] we shall describe two exact solutions of fj4~3l 
for rational e(z). In Section [7] we perform a perturbative analysis in the limit 
fi — ► of equation (|4"3"j) with e(z) given by fl4"2l . 

6 Toy models 

Equation (T4"3"l) can be solved exactly, if e(z) is rational. In this case U(z) 
is also rational, as follows immediately from (1451) . Furthermore, poles of 
U(z) in the complex z-plane can be located only at the poles of e(z). This 
information is sufficient to obtain the explicit solution of equation fj4~3]) . In 
the present Section we describe this procedure in more detail for two rational 



The reason for studying such exactly solvable "toy models" is twofold. First, 
on the exact solutions one can check the accuracy of perturbative methods. 
Second, function (HSl) . which arises in the Ising spin chain problem, can 
be approximated rather well by the function [B + Bie 2 (z)] with certain 
constants B , Bi, if h x is not close to the quantum phase transition point. 



In this subsection we consider a toy model, in which free fermions at h z = 
have the dispersion law ui{6) = — cos 9. Then function ffTTl) reduces to 




Kz + z- 1 ), 

(z + z- 1 ) - p(z 2 + z~ 2 ). 



6.1 Toy model 1 



e(0;6) = -2 cos(9/2) cosfl, 
and equation (|4ip takes the form: 



(52) 



[(-z - z- 1 ) cos(6/2) - E(Q)} d>(z) = - X z 




dz' <j){z') 
iri (z' — z) 2 



(53) 



After rescaling 



A 



E(Q) 



X 



(54) 



2cos(6/2) 



2cos(9/2) 
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we come to equation (pE31) with e(z) = e%(z). It is clear from (H5|) . that E/(z) 





a/ 




" 

T / 


\ z 















Figure 4: Integration paths in the z'-plane for the first toy model: arcs 71 
and 72 going from z = to z relate to (150]) ; loop C is the integration path 
in (155) . 

is analytical in the whole complex plane z for such e(z). Indeed, the left-hand 
side of ( I4"5j) is analytical inside the circle \z\ = 1, since g+(z) is analytical 
there, and e(z) is also analytical at \z\ < 1 except from the simple pole in 
the origin z — 0, where <7+(z) takes zero value. Similarly, the right-hand side 
of f|4"5]) is analytical outside the circle \z\ — 1. Therefore U(z) is a constant, 
and we can put U(z) = 1. The function is determined by (1501) . with 

2 = 0, C/(z) = 1, z ^(z) = -Alogz - -{z - z- 1 ). (55) 

We have set the initial integration point zq to zero in (1501) to provide g+{0) = 
0. Further, integration path in z'-plane should approach the point zq = 
from the left half-plane Imz' < to guarantee convergence of the integral 
in z' in the origin. Two integration paths 71 and 72, which satisfy all above 
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V 



Figure 5: Plot of Bessel function J„(l//x) versus v at fi = 0.5. 

conditions, are shown in Fig. HI Integration in (1501) along each path, either 
7i or 7 2 , should define t/ie same function g+(z). This requirement will be 
satisfied, if 



/(A) = 0, 



(56) 



where 



/(A) 



exp <{ - .F(z') 
Co 1 A* 



(57) 



and the integration path Co = 72 — 7i is shown in Fig. HI Condition (1561) 
determines the spectrum A for our first Toy model. After change of the 
integration variable z' = it -1 , the right-hand side in (1571) reduces (up to the 
sign) to the integral, which gives the well-known representation of the Bessel 
function, 



(0+) du 

exp[a(u — u )/2\ — 2iri J u {a), 

■ rvi If 



(58) 



with v = —X/n, a = Integration contour in (158|) "starts at infinity on 
the negative real w-axis, encircles the origin counter-clockwise, and returns 
to its starting point", see formula (5) in page 15 in [2]. So, the eigenvalues A„ 
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Figure 6: Dispersion curves E n (Q) for the first toy model for \ = 0-5 (left) 
and x — 0-2 (right) according to ( 1601) . 



of equation with e(z) = ei(z) are determined by zeros v n of the Bessel 
function J u (l/fi) as the function of its order v: 

\n = -v n fJ>, J Vn (l/}i) = 0. (59) 

A plot of J u (l/fi) versus u at fj, — 0.5 is shown in Fig. [5j For the energy 
spectrum E n (Q) in the first Toy model equations fl59|) . (jS"31) yield 

£„(e) = -^ x , (60) 

2cos(6/2)" 



where J Vn 



X 



0, n= 1,2,.... 



Fig. [6] shows the lowest dispersion curves E n (Q) for x = 0.5 and x = 0.2. 
Note, that E n {n) = n \ in this model. 

It is straightforward to extract the small-/x asymptotics from above exact 
results. In the fi —>■ +0 limit the integral in (|56|) is determined by the saddle 
points z± = —A ± iyl — A 2 of the function J-{z') . They lie on the circle 
\z'\ = 1, if |A| < 1. 

If |A| is far enough from 1, they provide two separate complex conjugate 
contributions to the integral in (15T|) . In the leading order in //, this yields 

A n o + sin0 o = 7T//( n - - ), cos0 a = -A„, (61) 



4 

in agreement with semiclassical quantization formula (T2 
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If A approaches —1, two saddle points z± merge at 1, and the integral 
(I57p becomes proportional to the Airy function [1] 

dOexpl - -(A+1)0+- j = 27r(2 / i) 1 / 3 Ai[-(A+l)(2) 1 / 3 /i _2/3 ] 

(62) 

in the limit fi —>■ 0. From (1561) we get in this case 

A„ = -1 + /i 2/3 2~ 1/3 ^ n , where Ai[-z n ] = 0, n = l,2, 

in agreement with (I3TI) . 

If A gets close to +1, z± approach —1. However, there are still two saddle 
point contributions to 1(E), coming from the the upper and lower edges of 
the contour Co at z' — — 1, see Fig. HI The resulting equation for A reads 

\nj Bi[(A-l) (2/^)1/3]' 

where Bi(x) is the second Airy functions, which solves (together with Ai(x)) 
the differential equation equation y"(x) — xy(x) = 0, and has the following 
integral representation [1] 

POO 

(3a)" 1/3 7r Bi[±(3a)" 1/3 x] = dt [exp(-at 3 ± xt) + sin(at 3 ± xt)]. (64) 

Jo 

Equation ( |63l) describes the crossover at A ~ 1 of the fi —>■ asymptotics 
from the semiclassical regime (|6TT) at A < 1 , to A n ~ fi n at A > 1 . 



6.2 Toy model 2 

In the first Toy model with ui(6) = — cos#, variation of the total momentum 
G only rescales e(9; G), see (|52l) . However, for the real free-fermion dispersion 
law ([9]), the topology of e(9; G) also changes with increasing G, compare Figs. 
[2] and [3j To mimic this property, let us modify the "free-fermion dispersion 
law" to 

u 2 (9) = -cos#-7cos(2#). (65) 
with positive 7. Then the function (fl71) takes the form 

e(0; G) = -2 cos(G/2) cos 9 - 2 7 cos(G) cos(2 0). (66) 
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Equation (l4"Tj) reads now 



{[-(z + z~ L ) cos(6/2) - (z 2 + z~ l ) 7 cos 9] - E(Q)]}<j>(z) (67) 
dz' </>(V 

-X z ' 



s vri (z'-z) 2 ' 



After rescaling 

g(g) X 7 cos 9 

cos(9/2)' ^ cos(9/2)' P cos(9/2)' 1 ' 

it reduces to (03]) with e(z) = e 2 (z) = -{z + a; -1 ) - p (z 2 + z~ 2 ). 

9ne can naturally come to this Toy model approximating the exact free- 
fermion dispersion law (Q by three initial terms of its Fourier expansion. 
This approximation is rather good, if h x is not too close to the quantum 
phase transition point h x — 1. 

It follows from (141)1) . that U(z) can be now written as 

U(z) = A X + A 2 (z + z- 1 ), 

with two constants Ai, A 2 . The function g + (z) determined by (150]) takes the 
form 



9+W 

where 



r dz' 










exp | 


Jo 







% -[T(z')-F{z)]}, (69) 



iT{z) = - X log z - ( z - - ) - £ ( z 2 - ^ 



2 / Z \ Z 2 



To provide convergence of the integral in (1691) , the variable z' should approach 
the origin inside one of two allowed sectors, which are defined by the condition 
Re(p/V 2 ) < 0, or explicitly 

7r/4 < argz' < 37r/4, and - 3tt/4 < argz' < -7r/4, if p > 0, 
— 7r/4 < argz' < 7r/4, and 57r/4 < arg z' < 37r/4, if p < 0. 

Definition ( 1691) becomes unambiguous and independent on the integration 
path between the points and z, if the constants A\ and A 2 solve two linear 
uniform equations: 

Wn(A) A 1 + W 12 (X)A 2 = 0, 
W 21 (X) A 1 + W 22 (X)A 2 = 0, 
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p >0 
(a) 



p<0 
(b) 



Figure 7: Integration contours C\ and C 2 in (1701) : 
Both Ci and C2 start and finish in the origin. 



(a) for p > 0, (b) p < 



where 



Wy(A) 



— - Wj(z') exp 



IV- 

0' 



1,2, 



(70) 



Integration contours Ci and C 2 for the two cases p > and p < are shown 
in Fig. [7J We fix the branch of T{z'} in (1701) by the condition = 0. 

Eigenvalues A n in the second Toy model are then determined by the re- 
quirement: 



W n {\ n )W 22 {\ n ) - W 12 {\ n )W 21 {\ n ) 



0. 



With known A n , we obtain finally E n (Q) by use of (!68|) . 
Note. In the same way one can find exact solutions of equation (|43|) with 
e ( z ) = J2k=i °k{ zk + z~ k ). In this case the eigenvalues A n of (H3|) are deter- 
mined from det W(X n ) = 0, where W (A) denotes a K x K-matrix. 
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The small-/! behavior of the second Toy model is more rich than in the 
previous case. It is determined by the location of four saddle points of the 
integrals in ([TO]) . These saddle points are the solutions of the equation e 2 (z) = 
A. We skip further discussions of the fj, — > asymptotics of the Toy model 
2, since it is very similar to that of original equation (1411) . which will be 
described in the following Section. 



7 Weak coupling expansion 



In this section we return to the original equation (HTj) . describing bound- 
spinons in the two-fermion approximation, and obtain its perturbative solu- 
tion in the weak coupling limit x ~~ * 0. Several perturbative schemes have 
been developed for the analogous problem in IFT [8, 10, 18]. Here we shall 
use a different procedure, which utilizes results of Section [51 
First, let us rewrite (jUJ) in notations of Section [5j 



[e(z) - X)(f>{z) = -nz 



dz' 



where 



e(z 



2\/h 



h, 



Si 



1 

ho- 



rn [z 



zv 



(71) 



-,1/2 



ZV 



exp(i6/2). 




h x 1 v, h x x v 1 , oo. Its 



/' 

Remind, that S\ is the unit circle, and <fr(l/z) = — 0( 
Function e(z) has six branching points 0, h x v , h x v ' 
Riemann surface C has four sheets £ Qj/ g with a,j3 = +, — , which are dis- 
tinguished by the signs of the first an the second terms in braces in (1721) at 
z = 1. To separate the sheets, we draw four cuts in the z-plane shown in 
Fig. [HI Equation (ITT]) is written in the sheet which we shall call "the 

physical sheet". 

Equations (Hlj) . fH51) define functions g±(z) and U(z) in £+,+. It follows 
from these equations and (172 p . that all singularities of U(z) in the physical 
sheet are the square root singularities of e(z). Therefore, U(z) is finite in 

£+,+• 

Let us turn now to the integral formula ([50]) . As in the previous Section, 
we should put z = in it to provide g + (0) = 0. The function T{z), de- 
termined by ([5T1) . is singular in the physical sheet at z — > 0, T{z) ~ z -1 / 2 . 



23 




Figure 8: Square root branching points of function e(z) determined by (1721) . 
and cuts on z-plane. 

Nevertheless, the integral in fl50|) converges, if the whole integration path lies 
in approaching to the origin either from the right, or from the left side. 

This follows from the fact, that e(z) is positive at real z in the physical sheet. 
Furthermore, three integral equalities must be satisfied 

/ — U{z) exp\-F(z)\ = 0, a = 0,1, 2, (74) 
Jc a z I A 1 J 

where the integration path Co is shown in Fig. HI and paths Cx,2 are drawn 
in Fig. [3(b). Note, that only two conditions in ( 1741) are independent, since 
Co — Ci + C 2 . Conditions (1741) are necessary to guarantee, that: 
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• the integral in (1501) defines a single valued function in the unit circle 

Si, 

• g+(z) given by (|5"0|) vanishes, when z approaches to the origin in £+,+ 
either from the left, or from the right side of the cuts, see Fig. [SI 

In the limit fi — > +0 the integrals in (1741) are determined by contributions of 
the saddle point of function J-(z), i.e. by the solutions of equation 



e(z 



A. 



(75) 



This equation has four solutions z a , Zf,, z a 1 , z b 1 in the Riemann surface C, 
which are given by 



Za Z 



A 2 cos(9/2) ± {[A 2 - 16sin 2 (G/2)][A 2 - 16 h\ sin 2 (6/2)]} 1 / 2 



8h x sin 2 (9/2) 



a 



a, b. 



Important for us are the saddle points, which lie on the circle \z\ = 1, or are 
close to it. In the simplest case, there are only two such solutions of (175]) : 
z a = exp(z6 l a ), and z~ l = exp(—i9 a ). Suppose, they are well separated, being 
far enough from 1 and -1. Then the saddle point expansion of the integral 
in (JT4"|) over the path Cq leads to the asymptotic expansion 



-F(Pa) = A*7r ( n — - ) + // arg (U(6 a + 56) exp 



P 



, (76) 



valid to all orders in /i. Here we use short notations T{9) = !F[z(9)], U(9) 
U[z(9)}, and 



AF{9 a + 89) = T{9 a + 



T[9 a ) 



59 2 . 



Averaging (f(S9)) implies formal term-by term integration 



(/(**)> 



1^(159/ (59) exp 


2 fj, 


f 00 ^ d59 exp 




2 /i 





25 



with f(89) expanded into a power series in 59. Under appropriate normal- 
ization of <fi(z), the function U(9 a + 59) can be expanded as 



U(9 a + 59) = l + J2J2 Ci ^ i5dl - 



(77) 



i=l 1=0 



To obtain the coefficients cu explicitly, one should write down the Neumann 
series, giving the formal solution of equation (ITT]) in the class of generalized 
functions, 



<j>{9) 
i<°>(0) 



(0) 



cf> {0 \9) + ^ 1 \9)+O( f /), 
C[5(9-9 a )-5(9 + 9 a )}, 



iir[e(9)-X] [sm 2 [(9-9 a )/2] sin 2 [(9 + 9 a )/2] 



(78) 



(79) 



Substitution of this expansion into (jUJ gives us a smooth function U{9) as 
a power series in ji 



U(9) = U {0) (9) + U {1) (9) + O(fi 2 ), 
U {0 \9) = ^[e(9)-e(9 a )\ 



(80) 



1 



Ll-exp[i(0-0«)] 1 - exp[i(0 + a )] 



which in turn can be expanded in 9 at 9 = 9 a . The result is brought to the 
form fl77J), if one puts C = -2m/e'(9 a ). 

The perturbation procedure described above should be modified, if equa- 
tion (1751) has four solutions z a , z b , z~ l , z^ 1 on the circle \z\ = 1, z a>b = 
exp[i9 a fi], see Fig. [31 If they are well separated, the zero-order function 
(f>(°>(9) in (1781 can be written as 



4> {0) (9) = C a [5(9 - 9 a ) - 5(9 + 9 a )\ + C b [5(9 - 9 b ) - 5(9 + 9 b 



(81) 



Using the leading saddle point approximation for the integrals in (17^1) with 
j = 1,2, we obtain a system of two equations 



u„ 



u a 



2lTfl 



-,1/2 



27T/X 



,1/2 



2lTfi 
27T/1 



lV2 



e m/<i+iF(e b )/n _ 
1/2 

e in/A-iT(6 b )/n _ g 
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for constants u a = U^(9 a ), = U^°'(6b), which are simply related with 
C a , Cb- Setting its determinant to zero, we get 

HOb) - ?{B*) = W\n ~ 5) + ( 82 ) 

in agreement with the semiclassical formula (129]) . 

It is easy to determine the eigenvalues A n of (ITTT) in the limit /i — > +0 for 
A > e{6 = 7r). Let us rewrite condition (fTIj) with a = as 

d0 Z7(0) exp[£F(0)//x] + / d0 Z7(0) exp[i.F(0)///] + 

-7T 7T 

/7T 
d0 Z7(0) exp[i\F(0)/fx] + 
-TV 

{l-exp[-iF(27r)//x]} / d6U{0) exp[iF{0)/n]. (83) 

J n 

In the second equality we have taken into account, that U(8 — 2tt) = U(8), 
and T{Q - 2tt) = T{Q) - T{2tx). 

If A > e{6 = 7r), the single saddle point z a in the integration path C is 
located in the interval (—1,0) on the real z-axis. Two contributions of this 
saddle point will cancel each other in fl83|) . if J~{2ti) = —27m fi with integer 
n, i.e. 



X n = — d6e(6)+iin. (84) 
^ Jo 

This asymptotic formula is valid to all orders^) in fj,. Rewriting it in original 
variables (1731 and using ffTTl) . one obtains the bound-spinon spectrum 

En{Q) = f d62u(6)+ X n (85) 
^ Jo 

for E n (Q) > e(ir, 9) and x -> +0- (86) 

This result has a clear physical interpretation. Condition flHBl) means, that 
a bound-spinon can be considered semiclassically as a large enough domain, 
bounded by two domain walls. These domain walls move back and forth 



'Corrections to (|84|) of order /i exp[— A(X)/fi\ can be obtained from 
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without collisions according to the classical equations of motion correspond- 
ing to the Hamiltonian (j!4p . The two terms on the right-hand side of (1851) 
give the kinetic and potential energy of domain walls averaged over their 
oscillations. So, the energy of such a bound-spinon does not depend on its 
momentum G, and hence, it can not move as a whole along the spin chain. 

The asymptotic formulas (ITS]) . (1521) . f[84"|) describe the small-/! behavior of 
eigenvalues A n of (1TT1) for generic values of A and G. Five crossover regimes 
are realized, when A approaches the critical values of e(6) and two or four 
solutions of equation (!75|) merge. All these regimes can be analyzed by means 
of almost the same perturbation procedure utilizing equations (1471) and (!74|) . 
Minor modifications in it are caused by further degeneracy of T{9) near the 
saddle point 9 : J-'iO) — ^{Oq) ~ {6 — 6> ) J where j = 3 or j = 5 instead of 
j = 2 in the case of Morse saddle point considered previously. In the [i —>■ 
limit this leads to perturbation expansions for A„ in powers of parameter /i 1 / 3 
or /i 1//5 . Below we present the leading terms of these expansions, skipping 
the details of calculations. 



1. 9 a ps 0, < G < Q m , where Q m = 2arccos/ia;, 
X n = e(0) + z n ^ 3 [e"(0)/2] 1/3 . 



(87) 



2. 6 a PS 7T, < < 7T, 

Ai 



tan 



Bi 



[A n -e(7r)] [-2/e"(7r)]V3^-V3 



[A n -e(7r)] [-2/6"(tt)]V3 m -2/3 



3. 6 m < 6 < 7r 

(a) 6 a ps # b 6*o > 0, where cos^o 



cos(e/2) 



cos(6 m /2) • 

A„,i = e(# ) + ^ 2/3 [e"(# )/2] 1/3 , 



A 



n.2 



where n — 1, 2, 
(b) b ps 0, # a > 0, 



cot an 



/i 4 



Ai 



[A n - c(0)] [-2/|e» 



Bi 



(89) 
(90) 



1/3^-2/3 



[A n -e(0)] [-2/|e"(0)|]V3 /U -2/3 



(91) 
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4. e = e m , e a « o, e b « z0 o , 



d 4 e(6)/d9 4 
6 



1/5 



C n > | e=o, 



(92) 



where e(0) = 2 (1 - Z^) 1 / 2 , 9 4 e(^)/^ 4 | e=0 = h*/ ^/T=W X , and n 
1,2 



The numbers z n and 2^ are consecutive zeroes of Ai(—z) and Ai'(-z) respec- 
tively, and c n are the solutions of equation 



dy 







sin y c n 

1 20 



exp 



x /x 5 
dx x 2 cos , 

20 



dyy" 



sin «c n 

1 20 



y 

h V C r 

20 y 



00 /x 5 
<ix cos 
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exp ( - — + y c n 



(93) 



ci = 1.787, c 2 = 3.544, c 3 = 5.086. 

Note, that all asymptotic regimes described above are realized also in the 
second Toy model, and the four regimes described by equations ( 1761) . ( 1541 . 
( 157|) . ( 15B1 are relevant to the first Toy model as well. 



8 Discussion 

In this paper we apply a technique developed in IFT to study the effect 
the discreteness of the lattice on the kink confinement in a non-critical one- 
dimensional system with explicitly broken Z 2 -symmetry. We study the model 
of the Ising spin-1/2 chain ferromagnet at zero temperature in the presence 
of a skew magnetic field, as a particular realization of such a system. The 
magnetic field h x , transverse to the easy z- axis, serves to induce the quantum 
phase transition at h x — ±1 and to allow kinks in the ordered phase \h x \ < 
1 to move. The longitudinal field h z breaks the Z2-symmetry and leads 
to confinement of kinks into pairs, the bound-spinons. We calculate the 
dispersion law E n (Q) of bound-spinons in the limit of small h z in the leading 
order in this parameter. 

The dispersion law of bound-spinons can be understood within a sim- 
ple heuristic picture, in which one treats a bound-spinon as bound state of 
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two classical particles (kinks) attracting each other with a linear potential 
proportional to h z . If these particles are close enough to each other, such 
that they can meat during their classical motion, then they can drift as a 
whole along the spin chain. In this case, one can determine the dispersion 
law E n (Q) of a bound-spinon from the Bohr-Sommerfeld quantization rule, 
if n > 1. 

On the other hand, large enough bound-spinons can be viewed as two 
independent well separated kinks. The linear potential leads to localization 
of an isolated kink in the discrete spin chain, similarly to localization of 
an electron moving in a periodic potential by a uniform electric field [24]. 
The classical motion of isolated kinks looks like oscillations around certain 
positions in the chain. If two kinks forming a bound-spinon do not meet 
during their classical oscillations, the velocity v n (Q) of the bound-spinon is 
zero. Since v n (Q) = dE n (Q)/dQ, we conclude E n (Q) is constant for such a 
bound-spinon. It is clear that the energy spectrum should be equidistant in 
this region 

E n+1 (Q)-E n (Q) = 2h z a, (94) 

since the bound-spinon state |$ n+ i(@)) contains one extra inverted spin in 
the domain bounded by two kinks, compared to the state |<3> n (6)). 

A more systematic theory is based on the perturbative analysis of the 
singular integral equation (jHJ), which is analogous to the Bethe-Salpeter 
equation in IFT . The equations have been obtained in the two-fermion ap- 
proximation and become exact in the limit h z — > 0. The bound-spinon 
energy spectrum E n (Q) is determined in this approximation as eigenvalues 
of equation (|4"T]) . We develop perturbation theory for this equation in small 
h z and, concentrating on its leading order, describe eight asymptotic regimes 
for the bound-spinon energy spectra E n (Q), instead of two regimes and 
flU) known in IFT. 

It was shown [8], that the IFT Bethe-Salpeter equation reproduces the 
energy spectra of stable "mesons" with reasonable accuracy not only in the 
limit h — > 0, but also at finite, and even at large values of the magnetic field 
h. If this situations holds in the discrete-lattice case, it would be reasonable 
to study the spectrum of equation fj4T|) at finite h z ~ 1 as well. Few steps 
in this direction give exact solutions of its simplified versions fl53|) and fl67|) . 
which were described in Section |6j 

In conclusion let us comment on the higher order corrections to the ob- 
tained bound-spinon spectrum E n (Q) in the limit h z — > 0. It is straightfor- 
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ward to determine the higher order terms in the weak coupling expansion of 
the eigenvalues of equation (14"T1) . as it was described in Section [0 However, 
one should also take into account corrections to equation (j4ip . which is an 
approximation by itself. First, the kernel G&(6, 8') in the integral in equation 
( |37|) has a regular part GQ Sg \6, 9') omitted in (|4T|) . Only few modifications in 
the perturbation schema of Section [7] are needed to apply it to equation (1371) 
with this term restored. Second, starting from the second order in h z one 
should take into multi-fermion corrections, which arise from the four-fermion, 
six-fermion, contributions to the bound-spinon eigenvector |$ n (0)). As 
we know from IFT [8,10,18], multi-fermion effects lead to renormalization of 
the fermion dispersion law uj{6) and coupling constant 2h z a, and should lead 
to the decay of unstable bound-spinons. All these effect are lost in the two- 
fermion approximation f[3"4"|) . but contribute to the bound-spinon dispersion 
law E n (Q) in higher orders in h z . A perturbative analysis of multi-fermion ef- 
fects in the discrete Ising model flTJ) is an interesting problem, which remains 
for future work. 
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